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The mechanism of incommensurate (IC) spin-density-wave (SDW) order observed in the 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase of CeCoIns is discussed on the basis of new 
mode-coupling scheme among IC-SDW order, two superconducting orders of FFLO with Big 
{d^2_y2) symmetry and 7r-pairing of odd-parity. Unhke the mode-couphng schemes proposed 
by Kenzelmann et al, Sciencexpress, 21 August (2008), that proposed in the present Letter can 
offer a simple explanation for why the IC-SDW order is observed only in FFLO phase and the 
IC wave vector is rather robust against the magnetic field. 
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Recent discovery of incommensurate (IC) spin- 
density-wave (SDW) order in the Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) phase of CeCoIns has attracted 
intense attention. -'^ This result is consistent with that 
of NMR Knight shift measurements suggesting the ex- 
istence of field-induced magnetism.^ Kenzelmann et al 
reported that an IC-SDW order in the c-direction with 
the wave vector Q = (hir/a, hir/a, 0.57r/c), h = 0.56 and 
a being the lattice constant in the a6-plane, exists in 
the so-called FFLO phase under the magnetic field H 
along the (1,1,0) direction. They also claimed that the 
wave vector Q is essentially constant against the mag- 
netic field strength H within an experimental resolution. 
In order to explain their results, they proposed a scenario 
on the basis of mode-couplings among the IC-SDW or- 
der Mq, d-wave pairing A^, additional superconducting 
orders with finite center-of-mass momentum with 
different combination of (i,n) on the symmetry require- 
ment. While the scenario explains the appearance of the 
IC-SDW order, it fails to explain the reason why the 
IC-SDW is observed only in the FFLO phase and the 
ordering wave vector Q is robust against the magnetic 
field. 

The purpose of this Letter is to present a simple and 
transparent explanation for these questions. It is cru- 
cial to examine an explicit form of the mode-coupling 
term on the semi-microscopic level beyond conventional 
group theoretical arguments which include ambiguities 
concerning combination of the wave vectors of order pa- 
rameters."'^'^ The mode-coupling term relevant to the 
present situation, in the GL region, is given by the Feyn- 
man diagram shown in Fig. 1, and the analytic expression 
is given as follows: 



Qo+qAL°Qo, 



(1) 



where Mq^+q, A^'^^ , and A^°q^ denote the IC-SDW mag- 
netization in the c-direction, the d-wave superconducting 
(SC) gap in the FFLO state, and an additional odd- 
parity "equal-spin" SC gap with finite center-of-mass 
momentum — Qo, Qo = (0.57r/a, O.dn/a, O.Stt/c), (the so- 
called TT-paring), respectively. The coefficient C in eq.(l) 



is given in terms of the Green function G'cr(k, ie„) of the 
quasiparticles in the normal state as follows: 

C = 9^9-Qo+ci/2TY,T.T. Mp)xliMP - Qo/2)xi°i X 

n p (T— ± 

Ga{p + q/2, ie„)Gff (p - Qo - q/2, ie„) x 

Gs(-p + q/2,-ie„), (2) 

where g's denote the pairing interactions for different 
center-of-mass momentum, and (j^eip) is the wave func- 
tion of d-wave pairing: e.g., 4>c{p) = (cos Pxa — cos Pya), 
and 0o(p) is the wave function of odd-parity state, xlr^ 
and denote the spin part of gap functions for even- 
and odd-party state, respectively. Here, the quantization 
axis of "spin" (label of crystalline-electric- field ground 
doublet) is taken as parallel to the magnetic field, i.e., 
in the (1,1,0) direction, and the direction of IC-SDW 
magnetization is along the (0,0,1) direction. Since xi-^ is 
antisymmetric with respect to interchange of a and a, 
i.e., XffCT = — xi-ff, and Xs's is not antisymmetric, the co- 
efficient C vanishes without the magnetic field: In other 
word, C is proportional to the magnetic field H or odd- 
function of i/. As a result, the mode-coupling term eq.(l) 
becomes the time-reversal invariant quantity. 

In order to explore the condition for the IC-SDW to be 
realized in practice, let us perform the analysis based on 
the GL free energy. For concise presentation, we denote 
A/q^+q as M and A^q^ as 4". Then, the GL free energy 
F is expressed as 



F = Fn 



4 



4 



(3) 
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where the coefficients um, a*, ^a/, and 6$ are assumed 
to be positive because both of M and ^ are assumed to 
be vanishing without the FFLO gap A^. The extremum 
conditions for M and ^ are 



dF 

= — ^ = qmM + buM^ 
dM 



(4) 
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Fig. 1. Feynman diagram for the mode-coupling term eq.(l) 
among IC-SDW magnetization, Mq^+q, d-wave SC gap of FFLO 
state, Aq, and the so-called vr-gap of odd-parity with center-of 
mass momentum Qq, A_qp. The solid line presents the Green 
function of the quasiparticles in the normal state. 
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Considering that CeCoIns is near the antiferromagnetic 
(AF) quantum critical point (QCP) with the ordering 
wave vector Qo without the magnetic field so that the 
IC-SDW whose ordering vector Q = Qo + q is also not 
far from the critical point, and anticipating the smallness 
of TT-gap ^, we neglect the cubic term of ^ in eq.(5) to 
obtain 

^ ^ -^M. (6) 

Substituting this into eq.(4), we obtain the equation for 
M: 



auM + bMM^ 



M = 0. 



This has M ^ solution if the condition 

is satisfied, giving the IC-SDW magnetization as 



om 



(CA, 



(7) 



(8) 



(9) 



The condition (8) will be satisfied if um is small 
enough and CAq is large enough. The former condition 
may be guaranteed by some circumstantial evidence that 
CeCoIus is located not far from the AF-QCP,**' '"^ and the 
latter condition is met in the FFLO state near the Hc2 
where the SC gap Aq remains to be large enough due to 
the first order transition at Hc2- The rather large value of 
Hc2 (comparable to the effective Fermi energy) supports 
the large value of the coefficient C. 

In the present model, the wave-vector of IC-SDW Q 
is not the same as that of FFLO q, but is given by 



Q = Qo + q 



(10) 



In the experiment of Ref. 1, the FFLO wave vector is con- 
sidered to be parallel to (1, 1, 0) so that the expected IC- 
SDW wave vector is along (0.5-|-q,0.5+q,0.5) direction in 
consistent with the experiment of ref. 1. However, if the 



FFLO vector was parallel to (1,0,0) or (0, 1,0), the ex- 
pected IC-SDW wave vector would be in (0.5 + g, 0.5, 0.5) 
or (0.5, 0.5 -|- q, 0.5) direction, respectively. This is a pre- 
diction of the present theory. 

If the magnitude of FFLO wave-vector q is given by 
q — 2^^H/hvp^ Q has some dependence on the magnetic 
field H. However, the relative effect is weaken by the 
presence of Qo in eq.(lO) in contrast to the case where Q 
itself is identified with FFLO wave-vector as discussed 
in ref. 1. Indeed, q is estimated in the FFLO phase as 
follows: 



2.4 X 10* (m"'). 



(11) 



where we have used iJ ~ 10 Tesla and wp — 7.5 x 10"^ 
m.^ Therefore, the ratio of q and Qo (= a/"") is 

q _ 4.6 X 10-1° 

Qo 71" 



X 2.4 X 10*^ ~ 3.3 X 10" 



(12) 



This value is about 3 times smaller than the observed 
value (0.56 — 0.5)/0.5 = 0.12. However, the Fermi velocity 
v-p used for the estimation in eq.(ll) is that near the 
SC transition temperature T^o at iJ = 0, and the mass 
enhancement at temperatures in FFLO phase is about 
5-6 times larger than that at T ~ Tco as estimated by 
the entropy balance argument. Therefore, it is reasonable 
that q, the deviation from the AF wave vector Qo, is 
given by the wave-vector of FFLO state. 

The robustness of IC-SDW wave vector Q against the 
magnetic field H is understood as follows: The range of 
variation of is 10.6 < i/( Tesla) < 11.4 so that the 
relative variation from the central value 7f = 11 is about 
3.6 X 10"^ This means the absolute value of variation 
of Q is 0.06 X 3.6 x 10"^ = 2.1 x 10^^ ^j^ich seems to 
be comparable to the error bar of IC-SDW wave vector 
although there is no error bar shown in the figure of Ref. 
1. Another factor, which reduces the effect of the mag- 
netic field H, is the _ff-dependence of the Fermi velocity 
which is an increasing function of H in general. 

The existence of the 7r-pairing of odd-parity can be de- 
tected, in principle, by the NQR measurement of the dy- 
namical susceptibility along the magnetic field (in (1,1,0) 
direction), which can give rise to an anomalous relax- 
ation of NQR spectrum as observed in Sr2Ru04^ with 
the same tetragonal symmetry as in CeCoIns. Indeed, 
since the wave function of 0o of the Cooper pair can be 
doubly degenerate in the odd-parity manifold, the spin- 
orbit coupling associated with relative motion of Cooper 
pairs can give rise to anomalous relaxation of NQR cor- 
responding to the oscillations of magnetization along the 
(1,1,0) direction.^ However, it may be left for future stud- 
ies to analyze in more detail how to detect the 7r-pairing. 

In conclusion, the mechanism of IC-SDW order in the 
FFLO phase of CeCoIns observed by neutron scattering 
has been clarified on the basis of a new mode-coupling 
scheme among IC-SDW order, two SC orders of FFLO 
with Big {d^2_y2) symmetry and 7r-pairing with odd- 
parity. The mode-coupling term proposed in the present 
Letter gives a simple explanation for why the IC-SDW 
order is observed only in FFLO phase and the IC wave 
vector is rather robust against the magnetic field. The 
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difference between the case of CeRhfns,® which is a sister 
compound of CeCoIns and IC-SDW state at H = 0, 
should be clarified by further investigations. 
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